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Abstract 

We consider a complex scalar field as SIMP dark matter in models with gauged Z 3 
discrete symmetry appearing as a remnant of dark local 17(1). Dark matter (DM) 
annihilates dominantly by the 3-to-2 scattering, due to the DM cubic coupling in 
combination with the DM quartic coupling or the Z' gauge and dark Higgs couplings. 
We show that a light Z' gauge boson makes DM in kinetic equilibrium with thermal 
plasma at freeze-out and it affects the DM relic density and perturbativity/unitarity 
constraints for DM self-interactions. We show that the large DM self-interactions are 
consistent with solving small-scale structure problems and explaining the DM halo 
separation recently observed in Abell 3827 cluster. Various bounds on the model from 
the SIMP conditions, DM self-interactions, Z' searches, DM direct/indirect detection 
experiments, and Higgs signals, are also discussed. 
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1 Introduction 


It is known that dark matter (DM) makes up about 85% of the matter density in the 
Universe [l] but the identity of dark matter has been regarded as one of the big mysteries 
in particle physics and cosmology communities. The Weakly Interacting Massive Particles 
(WIMP) paradigm has driven most of DM searches in terrestrial and satellite experiments, 
and typically, the abundance of WIMP is determined from the thermal number density 
at the time of freeze-out in terms of weak-scale DM mass and weak interactions with the 
SM particles. However, there has been no direct evidence for dark matter other than 
gravitational or cosmological effects and there are strong limits on the spin-independent 
DM-nucleon cross section |2,3|. Therefore, it is also important to look for alternative 
testable scenarios for dark matter other than WIMP. 


Strong Interacting Massive Particles (SIMP) |4',[5j have recently drawn attention due 
to the interesting property that a sub-GeV dark matter can be thermally produced due to 
(effective) 5-point self-interactions. SIMP dark matter might provide a new arena for the 
model building for dark matter that can be testable in the future dark matter searches. 
Large self-interactions of SIMP dark matter can play a role at the galaxy scales, for in¬ 
stance, for solving the cusp-core and too-big-to-fail problems |6, 7| and explaining the 
recently observed evidence for self-interacting dark matter in Abell 3827 cluster, namely, 
the separation between the dark matter halo and the galaxy’s stars (8 10 . Furthermore, 


the interactions of a light mediator between SIMP and the SM particles can be probed by 
various indirect and direct experiments [5]. 


In this article, as a SIMP dark matter candidate, we consider a complex scalar charged 
under a dark local I/(l)y, that is spontaneously broken down to a discrete Z 3 group, 
after another complex scalar gets a nonzero VEV. The remaining Z 3 symmetry ensures the 
stability of dark matter and generates the cubic self-interaction of dark matter dynamically 
due to the [/(l)y breaking. As a consequence, dark matter can self-annihilate through the 
3-to-2 scattering processes, thanks to the help of DM quartic coupling or Z' gauge and 
dark Higgs coupling. The dark Higgs could be a potential mediator to keep dark matter 
in thermal equilibrium via a mixing with the SM Higgs but it turns out to be insufficient 
due to the smallness of Higgs Yukawa couplings to the SM fermions at the time of SIMP 
freeze-out. Instead, the Z' gauge boson or dark photon can mediate for DM to scatter off 
the SM fermions sufficiently in the presence of a gauge kinetic mixing [5]. 


We show that the Z' gauge boson also gives sizable contributions to the DM annihilation 
and scattering processes for most of natural parameters of the model, while helping to 
realize SIMP as a thermal dark matter. We constrain the parameter space of our SIMP 


dark matter by considering the bounds coming from Bullet cluster 11 and halo shapes 12 


and find that our model can accommodate a potential hint for self-interacting dark matter 
from Abell 3827 cluster as estimated in Ref. 10 . Various limits from direct and indirect 


searches for a light Z' gauge boson are considered and compared to the SIMP conditions on 
the model. The direct detection of SIMP dark matter with electron recoil 13,14 and the 
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bound from indirect detection and the bound from Higgs signals are also briefly discussed. 

The paper is organized as follows. We begin to describe the model of a complex scalar 
dark matter with Z 3 gauge symmetry. Then, we present the results of the 3-to-2 an¬ 
nihilation and self-scattering processes of dark matter and consider the bounds on DM 
self-interactions. Next we show the complementarity of the SIMP conditions for the Z' 
searches at colliders and comment on the direct/indirect detection of dark matter and the 
impact on Higgs signals. Finally, conclusions are drawn. There are two appendices summa¬ 
rizing the Higgs-portal and Z'-portal interactions and the general formulas for annihilation 
and scattering cross sections. 


2 The model 


We consider dark matter as a complex scalar y having a charge q x = +1 under the local 
17(1) v symmetry, which is spontaneously broken by the VEV of another complex scalar (j) 
with charge = +3. Thus, the remaining discrete Z 3 symmetry is of the gauge symmetry 
and ensures the stability of scalar dark matter y. The SM particles are 


15,16 


origin 

assumed to be neutral under U(l)v- The U(l)y charges are summarized in Table 1. 



<t> 

X 

U(i)v 

+3 

+1 


Table 1: U(l)y charges. 


The Lagrangian for those scalars and the SM Higgs doublet H is given [] by 
C = -\v^ v - ~ sin iV^ + \D^\ 2 + \D, X \ 2 + \D,H\ 2 - V{</>,x,H) (1) 


where the held strength tensor for dark photon is V^ v = d^V u — d v V^ and covariant deriva¬ 
tives are = (<9 M - iq^goV^cj), D M y = - WxSP^y, and D^H = (d lt - ig'Y H B M - 

^igT a W^)H , and the gauge kinetic mixing between dark photon V/and hypercharge gauge 
boson B fl is introduced by sin £. Then, the dark photon communicates between dark matter 
and the SM particles through the gauge kinetic mixing. 


The scalar potential is V (</>, y, H) = Vdm + Vsm with 

H dm = -mJ|0| 2 + mJ|y| 2 + A^|0| 4 + A x |y| 4 + A 0x |0| 2 |y| 2 

+ ^0*y 3 + h.c.^ + A^|0| 2 |i7| 2 + A x fl-|y| 2 |i7| 2 , (2) 

Vsm = -m 2 H \H\ 2 + X H \H\ 4 . (3) 


1 The same local U(l)v model has been considered for WIMP dark matter in Ref. 
17(1) case has been also studied in the context of WIMP 


18 


or SIMP 


19 


17 


and the global 


dark matter. 
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Here, we note that A x , k are dimensionless self-couplings of dark matter, A^ x is the coupling 
between dark matter and dark Higgs, and A^#, X x h are Higgs-portal couplings. As will 
be shown in a later discussion on SIMP dark matter, the dark photon couplings to the 
SM fermions become important for keeping dark matter in kinetic equilibrium with the 
SM. Thus, although the Higgs-portal couplings could be interesting for direct/indirect 
detection, we focus on the Z'-portal in the later discussion. 

For a nonzero VEV of (j) with (0) = the U( l)y symmetry is broken to a discrete 
subgroup Z 3 and dark photon gets massive and can mix with the SM Z-boson. After 
expanding the dark Higgs as (j) = ~^{v' + h!) 1 the dark Higgs can mix with the SM Higgs 
by Higgs-portal interaction, X^h- The interaction terms for dark/SM Higgses and dark 
photon are given in the appendix A . 


3 Dynamics of SIMP dark matter 


When dark Higgs and dark photon are heavier than dark matter, a pair of dark matter 
annihilate into neither a single dark matter + a single dark Higgs (semi-annihilation) nor 
a pair of dark Higgses/dark photons (pair-annihilation). Nonetheless, dark Higgs and Z' 
contribute to DM annihilation and scattering processes as intermediate states. Dark Higgs 
and Z' can decay into a DM pair dominantly if kinematically allowed. For a sizable dark 
Higgs-DM coupling A^, x and Z' gauge coupling, dark Higgs and Z' would decay before the 
SIMP freeze-out, so they would not contribute to the observed DM relic density from their 
decays. 


Due to the remnant Z 3 discrete symmetry, dark matter self-annihilates dominantly 
through the following 3-to-2 self-annihilation processes, XXX* X*X* XXX XX *, and 
their complex conjugates, as shown in Figs, [l] and [2| The scattering processes can be 
mediated by dark/SM Higgses or dark photon/SM Z-boson. In the presence of Higgs 
and/or Z' portal interactions, the 2-to-2 annihilation such as XX* ff with / being 
the SM fermion is possible. The effect of the latter process will be discussed in the next 


section. We have adapted the CalcHEP package 20 to our model in order to compute the 


scattering matrix elements of DM annihilation processes. 

As will be discussed in a later section, the Z' portal interaction is important for keeping 
dark matter in kinetic equilibrium with the SM. We consider both Z' and dark Higgs 
interactions to DM in computing the 3-to-2 annihilation cross sections, but we ignore the 
Z' and Higgs portal interactions to the SM particles in this section. First, the squared 
amplitude for XXX* X*X* scattering is, in the non-relativistic limit, given by 

,o R 2 ( 9 200gpinl 

\M XXX * 12 - 11,702 D x 


>X X 


16m?, 

X. 


74A V - 117R 


m/ + m z , 


24 V m x( 3m x - 2m U “ A 0x( 43m x - 37m£ i )m|,/(9^) 


(4m£ 


+ m h J 


( 4 ) 
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Figure 1: Feynman diagrams for XXX* X*X*- 



Figure 2: Feynman diagrams for XXX XX* ■ 
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where R = \/2 kv' /(6m x ). Likewise, the squared amplitude for XXX XX* scattering is, 
in the non-relativistic limit, given by 


| M 


xxx-^xx* 


3 R 2 


rn: 


2A y + 9 R 2 + 


25 g 2 D m\ 


m £ + 


2A 0x m x( 13m x - 2 m 2 hi ) - Aj x (19 m\ - m 2 hi )m 2 z ,/(9g 2 D ] 


(9 rri* 


m 2 h Mml + m 2 ) 


(5) 


A nonzero cubic DM self-coupling only can make both 3-to-2 annihilation channels possible. 
Nonzero DM quartic coupling, dark gauge and dark Higgs couplings, can reduce or enhance 
the squared amplitudes, but the parameter space of DM self-couplings becomes strongly 
constrained due to unitarity bounds as will be discussed later. 

Assuming CP conservation in the dark sector, the Boltzmann equation for dark matter 
number density, tidm = n x + n x *, with n x = n x *, is 


dn 


DM 


dt 


+ 3L/udm 


= -4«' 
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av„ 


rel/XXX*—»X*X' 


+ Wei) 


XXX-*-XX' 


) (ri 


DM 


2 GO \ 
n DM n DMJ 


“Wei) 


xx*^ffy n DM 


in 


eq 

DM 


) 2 )- 


( 6 ) 


As a consequence, using the general formulas for the cross sections in Appendix B, the 
effective 3-to-2 annihilation cross section appearing in the above Boltzmann equation is 
obtained as 


(^rel) 


rel/3—>-2 


a (( Cry rel)xXX*->-X*X* + (^"Ael) 


xxx^-xx" 


) 


75 


15367 nnl 

X. 


| M 


xxx->xx 


+ I M 


XXX - ^xx* 


- “4 (7) 




Here, we note that the 2-to-2 annihilation channels in the dark sector are kinematically 
forbidden while mediators are assumed to couple to the SM particles weakly enough so 
that we can ignore the annihilation of dark matter into the SM particles. 

Therefore, when the 3-to-2 annihilation is the dominant process at freeze-out, the relic 
density is determined by n^ M (av 2 e ^^2 = H■> leading the relic density condition [4,[5] as 
follows, 

1/3 


m x oi e ff 


K 


Xi 


90 

g*n 


2 W 


( 8 ) 


where a e s is the effective coupling of dark matter in the 3-to-2 annihilation defined in 
eq. (Jt|) and T eq is the temperature at matter-radiation equality given by T eq = 0.8 eV. 
Taking xf — 20, and //* = 10.75 for the temperature 1 MeV < T < 100 MeV, and k = 
-frg* s c = 2.55 for c = 0.63g* : eq/.g*.s,e q = 0.54, the above condition becomes 


m x = 0.03 cieff (T^iHp) 1 ^ 3 . 


(9) 
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Figure 3: Feynman diagrams for XX* XX* and XX y XX ■ 

For a e ff = 1 — 10, we get m x = 35 — 350 MeV. In the next section, we search for the parame¬ 
ter space where a e g is of right magnitude for the SIMP dark matter, being compatible with 
perturbativity and unitarity and the bounds on the DM self-interactions. From the relic 
density condition (J9|) and the DM effective 3-to-2 interaction a e s given in eq. (J7|), we find 
one of the DM parameters dependent, i.e., the DM quartic coupling is now a dependent 
parameter by |A X | ~ oy^ 2 / R oc m x 2 / R. 


4 Bounds on SIMP dark matter 

In this section, we consider various bounds on SIMP dark matter, starting with bounds on 
DM self-interactions and Z'-portal couplings to the SM fermions. The SIMP conditions 
provide new constraints on the parameter space of dark photon, being complementary to 
the Z' searches at colliders. Direct detection of SIMP dark matter with electron recoil, 
indirect detection bounds and the impact on Higgs signals, are also discussed. 


4.1 Self-scattering cross section 


Dark matter self-scatters through the following 2-to-2 scattering processes, XX XXi its 
complex conjugate, and XX* XX *> as i n Fig. |3j Like in the 3-to-2 self-annihilation, we 
include both dark Higgs and Z' contributions to the scattering processes. 

First, the squared amplitude for the XX XX self-scattering is 


\M XX \ 2 


2 ^2A X + 3 R 2 + 


W D m 2 x 

m z> 


Kx m l ' \ 2 
$9 Hi 


( 10 ) 
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Figure 4: Left: Parameter space of m x vs R = \/2nv' /(6m x ), satisfying the relic density, 
the bounds from Bullet cluster and halo shapes. A^ x = 0.2, 0.4 is chosen for light to dark 
orange regions, respectively. We took mz> = 5 m x , go = 0.1 and m/ ll = 1.5m x . Right: 
Parameter space explaining the relic density and the separation between the dark matter 
halo and the stars of Abell 3827 with a se if/m x = 1 — 3 cm 2 /g [l0|. A^, x = 0, 0.1 is chosen for 
light to dark blue regions, respectively. In both cases, we took g D = 0.1, mz> = 5 m x and 
rrp M = 1.5m x , and perturbativity and unitarity are imposed. In both figures, we assumed 
a zero Higgs mixing angle, sin 0 = 0. Important constraints determining the remaining 
boundaries of the parameter space are written explicitly in the plots. 


On the other hand, the squared amplitude for the yy* —> yy* self-scattering is 

.2 


| M 


12 = 4 ( 2A - 9R 2 - 29d171x ' X *x( 2m x + m hi) m z' 


XX 


+ 


m|, 9^(4m 2 - m 2 hi )m 2 hi 


( 11 ) 


Therefore, in the non-relativistic limit for dark matter, the effective scattering cross section, 


— 1 ( -XX ,- X X i „XX \ .. -j-l, rXx _ ^xx 
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(\M XX ? + \M XX >\ 2 ). 


( 12 ) 


The self-interaction of dark matter is bounded from cr se if /m x < lcm 2 /g from Bullet cluster 


11 and halo shapes 12 


We also impose the perturbativity and unitarity bounds on the DM couplings, respec¬ 
tively, as follows, 


A x <4vr, \M xx \,\M xx A<^- 


(13) 
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Figure 5: Left: Parameter space of rriz' vs go, satisfying the relic density, the bounds from 
Bullet cluster and halo shapes. A^ x = 0.4, 0.6 is chosen for light to dark orange regions, 
respectively. We took R = 1.5, m x = 80MeV and m ^ = 120 MeV. Right: Parameter 
space of m/jj vs A^ x , satisfying the relic density, the bounds from Bullet cluster and halo 
shapes, go = 0.1,0.2 is chosen for light to dark orange regions, respectively. We took 
R = 1.5, m x = 80 MeV and rri z 1 = 400 MeV. In both figures, we assumed a zero Higgs 
mixing angle, sin# = 0. Note that the boundaries of the remaining parameter space are 
determined by unitarity in both figures. 


We note that XX ~t XX or its complex conjugate determines the unitarity bound domi¬ 
nantly. 


On the left panel of Fig. [4], we present the parameter space of DM cubic coupling k 
and DM mass by considering the relic density, the bounds from Bullet cluster and halo 
shapes as well as perturbativity and unitarity bounds. Here, we remind that the SIMP 
relic density condition, eq. §. was solved for A x , and the self-scattering cross section and 
perturbativity/unitarity bounds were written in terms of the remaining parameters. The 
bound on DM self-scattering leads to a lower bound on DM mass on the left of Fig. |4j On 
the other hand, the perturbativity and unitarity conditions cut off the parameter space 
with larger DM masses, because |A X | ~ alff / R oc m x 2 / R < An. We find that a relatively 
light dark Higgs with nonzero A^ x is necessary to satisfy the unitarity bound for cr se ii/m x < 
lcm 2 /g, because it cancels the DM self-coupling and dark gauge coupling contributions to 
the DM self-scattering amplitude. It turns out that the DM masses compatible with bounds 
on self-interactions are in the range between 20 MeV and 150 MeV depending on the dark 
Higgs and Z' couplings. On the other hand, the parameter space that explains the DM 
halo separation observed in Abell 3827 cluster with a se if/m x = 1 — 3cm 2 /g at 2cr 10 is 


also shown on the right panel of Fig. [4j In this case, the DM self-interaction required for 








<7 se if/m x < 1 cm2 /g, perturbativity <7 se | f /m x <1cm 2 /g, perturbativity 




m h i(MeV) /n h i(MeV) 


Figure 6: Parameter space of masses m/ ll vs mz>, satisfying the relic density, the bounds 
from Bullet cluster and halo shapes. We took m x = 80 GeV on left and 120 MeV on 
right. In both figures, gr> = 0.1, 0.2 is chosen for light to dark orange regions, respectively, 
and we took R = 1.5 and A^ x = 0.4. We assumed a zero Higgs mixing angle, sin 0 = 0. 
Important constraints determining the boundaries of the remaining parameter space are 
written explicitly in the plots. 
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Abell 3827 can be accommodated for A^ x = 0, namely, without dark Higgs contribution, 
and the DM masses range between 30 MeV and 60 MeV. 

Keeping the relic density, the bounds from Bullet cluster and halo shapes as well as 
perturbativity and unitarity bounds, we have searched for the parameter space in dark Z' 
mass vs gauge coupling on left and dark Higgs mass vs \^ x on right, in Fig. [5] We find 
that Z' mass is bounded from above for a given Z' gauge coupling while a relatively light 
dark Higgs with sizable X 0x is preferred. Lastly, in Fig. [6j we also depict the correlation 
between Z' and dark Higgs masses, depending on Z' gauge coupling and DM mass. In 
all Figs. 4-6, we marked the important constraints determining the boundaries of the 
remaining parameter space in the plots. 


4.2 Z' portal couplings to the SM 


We need a nonzero interaction between dark matter and the SM particles in order for dark 
matter to keep in kinetic equilibrium until freeze-out, while the consequent DM annihilation 
into a pair of the SM particles should be subdominant as compared to the 3-to-2 DM 
annihilation. These are the SIMP conditions. The DM interactions with the SM particles 
can be induced by Higgs or Z' portals, leading to DM annihilation or scattering as in 
Fig-0 In the SIMP scenario, for a sub-GeV DM mass, it was shown that the Higgs- 
portal coupling is not sufficient for keeping DM in thermal equilibrium j5,19 . But, for 
clarification, we present a separate discussion on the Higgs portals in a later section. In 
this section, we focus the Z'-portal couplings and discuss bounds from SIMP conditions 
and Z l direct searches. 

In the presence of gauge kinetic mixing, dark matter can self-annihilate into a pair of 
the SM fermions through dark photon exchange in the s-channel as shown in the third 
Feynman diagram in Fig. 7. The amplitude for X.(Pi)x* (P 2 ) —>■ f{ki)f(k 2 ) process is 


LM ann = ~ig D (pi ~ P2T 


-IQ [iv 


s — mi + imz'T z> 


-i)eQ f e u(k 1 )Yv(k 2 ) 


(14) 


where Qf is the electromagnetic charge of the SM fermion and e ~ cos 6w£ f° r m z> mz- 
Then, the squared amplitude is given by 

2 4 e Q y& 

K vl ann| / 2 \2 

(s — m x y 

Therefore, in the non-relativistic limit for dark matter, the annihilation cross section for 
XX* ff with Q f = -1 for charged leptons is 


^(ki-{pi-p2)){k2-{pi-P2))-(ki-k 2 ){p 1 -p2) 2 ~m 2 f (pi-p2) 2 ■ (15) 


(o'^rel)c 


dn(^l\M 


647rp'j>>2 


2 £ 2 ‘e 2 g 2 D m 


7r[(4m£ — m 2 z ,) 2 + \m x J m 


T\_5 2 


(16) 
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Figure 7: Dark matter annihilation and scattering through Higgs-portal (Left panel) and 
Z'-portal (Right panel). 


On the other hand, dark matter can scatter with the SM fermions that are in thermal 
equilibrium at dark matter freeze-out. The amplitude for X.{p\) f {p'l) —> x(ki)f{k 2 ) process 
as shown in the last Feynman diagram in Fig. 7 is 


iM scM = ~ig D (Pi + hY 




t — rn 2 z , 


(- i)eQ f eu(k 2 )YYP2 ) 


(17) 


with t = (pi — ki) 2 . Then, the spin-averaged amplitude squared is given by 
2j £^r 

|7W SC att| J = — - f 2 2(k- 2 -(pi+ki))(p- 2 -(pi+ki)) — (k 2 -p 2 )(pi+ki) 2 +mj(pi+ki) 2 . (18) 

[t m z ,) l - J 

Around the freeze-out, the SM fermions in thermal equilibrium are relativistic so we can 
take the momenta relevant for the DM-SM fermion scattering to be p\ = (m x ,p,0),p 2 = 
(p, — p, 0) and k\ = ( m x ,pcosd,psmd ), k 2 = (p, — pcosd, — psind) in the center of mass 
frame where 9 is the scattering angle of dark matter. Therefore, we obtain the kinetic 
scattering cross section for xf —> xf with Qf = — 1 as 


(cr^ei) 


scatt 


64:7Tpip2 J 

3 e 2 e 2 g 2 D m 


dn^\M scatt \ 2 ) 


271171 %, 


m. 


rri 2 


(19) 


The cross section for xf xf process is given by the same formula in eq. (19). 
The SIMP conditions are imposed at freeze-out temperature, as follows, 
^DM(c’’u re i) ann < nl M (av 2 el ) 3 ^ 2 < nsM(c r n r el)kin- 


( 20 ) 


The above conditions (20) constrain the 2-to-2 annihilation cross section and the dark 


matter kinetic scattering in terms of 8i t2 defined in eqs. (16) and (19) and a e s, as follows, 

Si < 2.4 x IQ’ 6 a e s, 8 2 > 10 ~ 9 a^. (21) 


Then, using eqs. (16) and (19) with the SIMP conditions bound the dark gauge 


coupling and the gauge kinetic mixing as a function of dark photon mass as 


2.2 x 10- 8 «;/ 2 ^ < \e\g D < 4.4 x 10- 5 « eff ^ 


mi 


mi 


+ n, 


mi 


z 'm 2 
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( 22 ) 
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m x =60MeV, go=0.3 m x =60MeV, g D - 0.6 




Figure 8: Various limits on the parameter space of gauge kinetic mixing £ = £ cos 9w vs 
mz'- We have chosen m x = 60MeV and gr> = 0.3(0.6) on left(right). SIMP I corresponds 
to the region that does not fulfill the condition, tidm {crv Te \) an n < n T>M{ aV T&)z-*2i while 
SIMP II corresponds to the region that does not satisfy the thermal equilibrium condition, 
n DM (^1)3^2 < n SM (crU r el)kin- 


The parameter space of dark photon can be constrained by direct searches at colliders 


21-27 and lepton g — 2 experiments |28J, etc. There has been a recent study on the limits 
on light dark matter in the context of a vector messenger 27 . In Fig. [8j various limits 


are imposed in the parameter space of the gauge kinetic mixing e ~ £ cos 6w vs mz', and 
some values of DM-electron scattering cross sections are presented. The direct limits from 


monophoton + MET at BaBar in light blue 22-25 and the future prospects at Belle2 in 
dashed lines 


24 


are given; the indirect limits from electron g — 2 in orange and EWPT 
are presented; the muon g — 2 favored region in green is shown too. The SIMP conditions, 
(SIMP I in red is the dominance of the 2-to-2 annihilation into the SM fermions while 
SIMP II in purple is the decoupling of SIMP from thermal plasma), rule out the parameter 
space that is not covered by direct and indirect limits on dark photon. The parameter 
space with < m x in gray is ruled out due to the dominance of the 2-to-2 annihilation 
into dark photons. 

The recent BaBar search for monophoton + dilepton (not shown) 126] can limit e at 
the level of 10 -4 — 10 -3 for dark photon masses in the range 0.02 — 10.2 GeV, but the 
limits depend on the decay branching fraction of dark photon. When the decay of dark 
photon into a pair of dark matter is open, it dominates easily, so the previously mentioned 
monophoton + MET is more constraining. For instance, the monophoton + MET at 
BaBar applies for dark photon mass being above 120 MeV in Fig. [8j while the monophoton 
+ dilepton at Babar applies and limits e at the level of 10 -3 for 20 MeV < mz> < 120 MeV. 
There are other bounds from beam dump experiments (not shown) j 21,25 that limit e at 
the level of 10 -3 in our parameter space below mz' — 0.1 GeV. The planned SPS target 
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Figure 9: Direct detection cross section between dark matter and electron as a function of 
Z' mass. We have chosen £ = 1CT 4 , and go = 0.1, 0.3, 0.6 for blue dashed, orange dotted 
and purple dot-dashed lines, respectively. The cross section depends little on m x in the 
parameter space satisfying the SIMP relic density. 


experiment at CERN 29 could improve the limits from beam dump experiments further. 


4.3 Direct detection 


In the case of m e ,m x ,mz' p — m x v dm, which is true of our SIMP dark matter [^J we 
can take the momenta, p\ = (m x ,p, 0),p 2 = (m e , — p, 0) and k\ = (m x ,pcos9,psm9), 


k >2 = (m e , — pcosO, — psind) in the calculation of the matrix element squared in eq. (18). 


Then, the DM-electron scattering cross section with Z'-portal interaction is, at present, 


1 


odd = 


647t 2 s _ 
£ 2 e 2 g 2 D p 2 
irm%, 


dtt jhj |A1 sca tt| 


(23) 


where p, = m e m x j (■ m e +m x ) is the reduced mass of the DM-electron system. For m x m e , 
the DM-electron scattering cross section depends little on the DM mass, but is sensitive 
to the gauge coupling and mass of dark photon. The electron recoil energy can be used to 


search for the signal of a sub-GeV light dark matter 113,14 


2 When the momentum transfer of order m x v dm is greater than Z' mass, then the DM scattering cross 
section becomes enhanced due to a Coulomb singularity 14 . 


13 












In Fig. [9j we depict the DM-electron scattering cross section as a f un ction of the Z' 
mass by varying the dark photon gauge coupling g D between 0.1, 0.3, 0.6 but with fixed 
gauge kinetic mixing £ = 10 -4 . The result is insensitive to the SIMP mass, which ranges 
between 30MeV and 150 MeV as shown from the left figure of Fig. |4j In Fig. [8j we also 
present the DM-electron scattering cross section on the parameter space of e vs mz>- The 
direct search for dark matter with electron recoil in XENONIO [13] can give stronger 
bounds than beam dump experiments, at the level of e rv_/ 10 1 (10 3 ) for rriz' ~ 0.06 GeV 
with gp, = 0.3(0.6). XENONIO 13 constrains the DM-electron scattering cross section 


to be below 2 x 10~ 36 cm 2 at best at about 30 MeV so the SIMP parameter space in our 
model that is compatible with direct Z' searches satisfies the current XENONIO limits. 
But, the case with a light Z' gauge boson could be also probed by future superconducting 
detectors 


14 


4.4 Higgs portal couplings and Higgs signals 


We discuss the impact of Higgs portal couplings on the DM scattering and annihilation 
processes. As shown in the second Feynman diagrams in Fig. 7, the SM Higgs and the dark 
Higgs contribute to the DM scattering with the SM charged leptons. The corresponding 
annihilation cross section is given by 




1 { m fV 

167rm?, \ v J 

X 


Vhix*\ 


+ 


Vh 2X*X 



(3 m x Tp + mj) 


(24) 


where Tp = m x / 20 is the freeze-out temperature and yh lX *x = s i n $(A ' cos 9 — \ x hv sin 9) 
and Vh 2X * x = — cos 6(\<t> x v' sin 9 + A x hv cos 9). 

On the other hand, the annihilation cross section for XX* f f is 


{<7V r el) xx *^yff ~ ^ 




m 2 f \ 3/2 

ml / 


Vhi 


xx 


4m2 


+ 


Dh 2 X *x 
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hi 


4 


m 


h 2 


(25) 


We digress to the bounds from Higgs signals, which constrain the Higgs mixing angle 
and the DM coupling to the SM Higgs, X x h- The bound on the branching fraction of the 
Higgs invisible decay is BR(/z 2 —> XX*) < 0-29 at 95% CL from the VBF Higgs production 
at ATLAS 1331. The ZH production at ATLAS |34j and the VBF +ZH production CMS 
leads to less strong bounds, BR(/r 2 —> XX*) < 0-75 and BR(/r 2 —>■ XX*) < 0-58, 


35 


respectively. On the other hand, the Higgs signal strength is bounded to /r > 0.81 at 95% 
CL from ATLAS/CMS data combined (36). 


The partial decay rates for additional Higgs decay modes in our model are 


r(/i 2 ^ xx*) = 


9h 2X * X 

167 rm /,. 2 


4- 


4m^ 


m 


(26) 


h2 
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KhV 2 

327rmh 2 
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Am 2 hi 


771 


(27) 


h2 


Then, from the bound on the Higgs invisible decay, we obtain \Vh 2 x*x\/ v ~ 0.011. On 
the other hand, scaling the total and individual Higgs decay rates as well as the Higgs 
production cross section approximately by cos 2 9 , the Higgs signals bounds the Higgs mixing 
angle as sin 9 < 0.44. In our case, the Higgs mixing angle is approximated to sin 9 ~ 
X^hvv'/ rn 2 hrt < 0.016 for rri} l2 m^, v' ~ 1 GeV and X^h < 1. Thus, it is consistent 
with the Higgs visible decays. Given the bounds, \yh 2 x*x\/ v ~ 0.011 and sin 9 < 0.016, we 
find that the dark Higgs contributes dominantly to the DM scattering and annihilation, as 
compared to the SM Higgs. As for Z' portal, we can parametrize the 2-to-2 cross sections 


with dark Higgs exchanges by (crv) ann = h 2 /m 2 and (cry) scatt = <5f/m 2 , in eqs. (25) and 
(24), respectively. Therefore, taking the charged lepton to be electron or positronjj and 
choosing A^ x ~ 1, we get hi ~ <5 2 ~ 10~ 10 . Therefore, unless A^ x > 10, the Higgs portal 
does not satisfy the SIMP condition. 


4.5 Indirect detection 


The SIMP dark matter can be searched for by indirect detection experiments with diffuse 
gamma-rays. Although the Fermi-LAT is sensitive to gamma-rays of energies from 20 MeV 
to > 300 GeV, its current limits, for instance, from dwarf spheroidal galaxies (dSphs) of 


the Milky Way 30 , do not extend down to 20 MeV. On the other hand, the AMS positron 


fraction in primary cosmic rays is within the energy range from 0.5 to 500 GeV [31]. The 
current best upper limits on e + e _ {/i + y~) annihilation cross section is 5 x 10 _27 (10 _26 )cm 3 /s 
at m x = 5 GeV from Fermi-LAT dSphs 30 and 10~ 28 (10~ 27 )cm 3 /s at similar DM masses 


from AMS-02 positron fraction 32 . More importantly, there are diffuse X-ray and gamma- 


ray observations near the center of the galaxy such as INTEGRAL [371, COMPTEL 38 and 


EGRET 39 and Fermi-LAT 40 . In particular, COMPTEL leads to the bounds stronger 


than the thermal annihilation cross section around 3 x 10 26 cm 3 /s, for DM masses below 


100 MeV 41 


SIMP can annihilate into e + e or h + /j, depending on its mass, through Z' or Higgs 
portal with the annihilation cross section, 


Wv rel )i+i- = (l.2 x 10 27 cm 3 /s) ^ 


hi \ 2 /100 MeV x 2 


10 


-6 


m , 


(28) 


In the case of Z' portal, the annihilation cross section is p-wave suppressed as hi(T) = 
5^fx/tJTf where = hi(Tp) < 10” 6 from the SIMP condition, we get hi(T) = (0.01 — 


3 When DM is heavier than muon, it can annihilate into a muon pair too. But, the DM scattering rate 

with muon is suppressed by Boltzmann factor at freeze-out [5], because the typical freeze-out temperature 

is small, Tp ~ m x /20 < 20 MeV. 
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0.1)(5i ; f 10 -8 — 10 - ' at present. On the other hand, the annihilation cross section with 

Higgs-portal is given by 8\ = 10 _10 (m;/m e ), which could be as large as 10” 8 for mi = m^. 
Therefore, in either Z' or Higgs portal, there is no current bound from indirect detection. 


5 Conclusions 

We have shown that models with gauged Z 3 symmetry provide a consistent SIMP scalar 
dark matter with sub-GeV masses that can be in kinetic equilibrium with the SM fermions 
via Z'-portal interactions. The resulting large self-scattering cross section of the SIMP 
dark matter could explain the DM halo separation observed in Abell 3827 cluster, although 
there is a tension with the previous bounds coming from Bullet cluster and halo shapes. 
Inclusion of Z' gauge boson and dark Higgs changes DM annihilation and scattering pro¬ 
cesses and opens up a new parameter space of the SIMP self-couplings and masses. We 
found that SIMP conditions are complementary to indirect and collider Z' searches, DM 
direct/indirect detection and Higgs signals. The SIMP parameter space could be explored 
further by future collider and indirect/direct detection experiments. 
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Appendix A: Higgs and 
matter 


Z' portal interactions to dark 


After the breakdown of clectroweak symmetry and dark U{l)y, we can take the SM 
Higgs and dark Higgs in unitary gauge, 

On the other hand, the VEV of scalar dark matter \ is assumed to vanish. Then, the mass 
matrix for h', h is given by 


f 2\ (j> v' 2 A qhv'v \ 

V A qhv'v 2\ h v 2 ) ' 


(A.2) 
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where the VEVs are determined by the following conditions, 

2 _ 4A 4 m 2 H - 2A^m| /2 _ 4A- 2\ (j)H m 2 H 

4A^A^ - A| h ’ 4A//A0 - A^ H 

Then, the conditions for a local minimum are 

4A^A# > A 2 H , A^, >0, A h > 0, 

2X,pm 2 H - A> 0, 

2A H m\ - \<pH j m 2 H > 0. 

The mass matrix can be diagonalized by the rotation by 

/ hi \ _ f cos $ — sin 9 \ / h' A 
V ^2 / V cosd y \ h y 

where hi, h 2 are mass eigenstates. The mass eigenvalues of Higgs-like states are 


m 


2 

h\ ,h ,2 


and the mixing angle is 


A ^' 2 + X H v 2 =F \J (A 0 u /2 - X H v 2 ) 2 + X 2 0 H v , 2 v 2 


tan 29 


X^hv'v 

X H v 2 - X^v' 2 ' 


(A.3) 


(A.4) 

(A.5) 

(A. 6 ) 


(A. 7) 


(A. 8 ) 


(A.9) 


Consequently, the self-interactions and Higgs-portal interactions of dark matter are 


^scalar = {Xfav' cos 9 - X xH v sin 9)hi \x | 2 + (A^v' sin 9 + X xH vcos 9)h 2 \x \ 2 

+ ^(A</>x cos 2 9 + A xff sin 2 6>)/ii|x| 2 + ^(A 0X sin 2 6» + A XjF / cos 2 6»)/i||x| 2 

+ (A 0 X - A x #) sin6*cosd/?ih 2 |xi 2 + ^/tcosd hiy 3 + +^Ksin 9 h 2 x 3 + h.c.j 


+ A xlx| 4 + (gj«u , X 3 + h.c.). 


(A. 10 ) 


On the other hand, the gauge kinetic and mass terms can be diagonalized by 


Bfj, \ / cw —sw c ( + ~ — ^ c C \ / A /t A 

bh 3 = % cws<; I Z\ t , (A.11) 

V n ) \ 0 -s c /c c c f /c ? / \ z 2jl ) 


where (. B kh 3 , V^) are hypercharge, neutral-weak and dark gauge fields, and (A M , Zi M , Z 2m ) 
are mass eigenstates, and s w = sin 6 ^, cw = cos 6 ^, etc. The mass eigenvalues of Z-boson 
and dark photon are 


2 x 
m, , = - 

b 2 2 


m 2 z ( 1 + Swt}) + my/q ± J (m|(l + s 2 w tf) + m^/c |) 2 - 4m|m^/c| (A. 12 ) 
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where m z = |( g 2 + g' 2 )v 2 and rriy = 9g^v' 2 , and the mixing angle between Z-boson and 
dark photon is given by 


tan 2^ = 


m 2 z sw sin 2£ 


iriy — m 2 z ( 


^2 — o2 o2\ 


(A-13) 


As a consequence of the basis rotation, the current interactions can be written as 


£ = 


e A^EM 


+-Zp t 
-\-Z 2 n 


eeJ em + 2 S wc w ^ ~ tw£ / t ^> J z ~ 9d ^ J d 

— eeJ^ M + — (s^ + t\ye)J z A 9d~Jd 

ZS\yCw Cg 


(A.14) 
(A.15) 


where £ = and J^and J/, are electromagnetic, neutral and dark currents, 

respectively. For my <C m^, we get £ ~ -s^, so the neutral current interaction of dark 
photon is negligible due to s ^ + t\y£ — C + — 0. 


Appendix B: General formulas for cross sections 


The cross sections for 3-to-2 self-scattering, DM-DM annihilation, DM-SM kinetic scat¬ 
tering and 2-to-2 self-scattering, in order, can be written in terms of the squared amplitudes 
as follows, 


where 

states 


(av 2 cl ) 3 ^ 2 

((TWrel)ann 
(<™rel) scatt 


75 


3847t?t?T 

X, 


| M 


3->2 


^"self — 


327T m 2 
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167 im? 
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647 Tin 2 
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| M 
| M 


2 

ann | i 


scatt | j 


| M 


self 


(B.l) 

(B.2) 

(B.3) 

(B.4) 


a 1/2 factor is taken into account in the squared matrix elements for identical final 
and we have ignored the SM particle masses. 
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